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Abstract. We propose to use null vectors in conformal field theories to derive 
model Hamiltonians of quantum spin chains and corresponding ground state wave 
function(s). The approach is quite general, and we illustrate it by constructing a 
family of Hamiltonians whose ground states are the chiral correlators of the SU{2)i^ 
WZW model for integer values of the level k. The simplest example corresponds to 
k — \ and is essentially a nonuniform generalization of the Haldane-Shastry model with 
long-range exchange couplings. At level k = 2,we analyze the model for TV spin 1 fields. 
We find that the Renyi entropy and the two-point spin correlator show, respectively, 
logarithmic growth and algebraic decay. Furthermore, we use the null vectors to derive 
a set of algebraic, linear equations relating spin correlators within each model. At level 
k = 1, these equations allow us to compute the two-point spin correlators analytically 
for the finite chain uniform Haldane-Shastry model and to obtain numerical results for 
the nonuniform case and for higher-point spin correlators in a very simple way and 
without resorting to Monte Carlo techniques. 
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1. Introduction 

Chains of interacting spins are important models for investigating properties of 
quantum many-body systems. The fact that the dimension of the Hilbert space grows 
exponentially with the length of the chains does, however, limit the applicability of brute 
force numerical computations. This makes analytical insights particularly valuable, 
both as a direct source of information and as a starting point for constructing and 
testing approximate numerical methods. Important analytical techniques include the 
Bethe ansatz [1] and the construction lying behind the AKLT model [2]. The former 
was introduced by Bethe in 1931 to solve the antiferro magnetic Heisenberg model with 
nearest-neighbour interactions and has turned out to be a more generally applicable tool 
to solve particular spin models, and the latter has led to a number of models for which 
both the Hamiltonian and the corresponding ground state(s) are known. Solvable models 
with long-range interactions are more rare, but the Haldane-Shastry (HS) model [3l H] 
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is an example. In this model, spin 1/2 particles are positioned uniformly on a circle, 
and each spin interacts with all other spins through a two-body exchange interaction 
whose strength is inversely proportional to the square of the distance between the spins 
measured along the chord. 

In the present paper, we propose a new method to find matching Hamiltonians 
and ground state wave functions, which naturally leads to models with long-range 
interactions. The approach, which is presented in detail in section |2l relies on the 
existence of null vectors in conformal field theories (CFT), and we shall specifically 
consider the SU{2)k WZW models, where k is a positive integer known as the level of 
the algebra. The vacuum expectation value of a product of primary chiral conformal 
fields is zero if one of the fields is a null vector, and we use this fact and the Ward 
identity to derive equations of the form Hip = 0, where if is a Hermitian and positive 
semidefinite operator and ip is a chiral correlator (i.e., the vacuum expectation value 
of a product of primary chiral fields). Explicit analytical expressions for several of the 
relevant chiral correlators have been found recently in [5]. 

In addition to deriving new spin models, the null vectors enable us to derive a set of 
simple algebraic linear equations, which relate different n-point spin correlators within 
each model (i.e., expectation values of the form {ip\t°l^ . . .t^"\ip) , where t"^^. is the a*'^ 
component of the spin operator acting on spin number ij). The approach turns out to 
be particularly valuable for /c = 1, where it leads to a closed set of equations for the two- 
point spin correlators, a closed set of equations for the two-point and the four-point spin 
correlators and so on. As detailed below, this allows us to obtain a number of results 
for spin correlators either analytically or through a very simple numerical computation. 

After explaining the general approach, which is valid for all k, we turn to a more 
detailed investigation of models obtained at levels k = 1 and = 2 in sections |3] and 
m respectively. In the SU{2)i WZW model, there are only spin and spin 1/2 chiral 
conformal fields, and since the fusion properties of the spin field are trivial, we find the 
explicit expression for a Hamiltonian which has the chiral correlator of A^ spin 1/2 fields 
as its ground state in section 13.11 The Hamiltonian involves only two-body interactions, 
and in section 13. 2[ we show that it is essentially a nonuniform generalization of the 
HS model with long-range exchange interactions. We also compute the spectrum for a 
number of choices of distributions for A^ = 6 and discuss symmetries and the origin of 
degeneracies. The Hamiltonian of the nonuniform HS model was put forward recently 
in [6], where it was shown analytically that the chiral correlator of A^ spin 1/2 fields is 
an eigenstate of it. This was done by using the Knizhnik-Zamolodchikov (KZ) equation 
[7], and we provide the details of that derivation in Appendix A In contrast to the 
present work, however, this approach cannot be straightforwardly generalized to higher 
values of k, and it does not automatically guarantee that the chiral correlator is, indeed, 
the ground state. 

The nonuniform HS model was compared to other spin models through exact 
numerical computations for small A^ in [6], and the finite chain Renyi entropy and 
the two-point spin correlation function was studied numerically through Monte Carlo 
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simulations. In the present paper, we instead concentrate on the analytical and semi- 
analytical expressions for the spin correlators that follow from the properties of the 
null vectors through the aforementioned linear equations. We rederive the expression 
for the two-point spin correlator of the uniform HS model in the thermodynamic limit 
and furthermore demonstrate that the two-point spin correlator can also be computed 
analytically for the finite chain uniform HS model. These results are presented in section 
13. 3[ where we also provide further numerical results for two-point and four-point spin 
correlators. 

Section m is devoted to an investigation of the case k = 2. In this model, there are 
spin 0, spin 1/2 and spin 1 chiral conformal fields, and we shall, in particular, analyze 
the spin chain model in which the wave function is the chiral correlator of a product 
of spin 1 fields. In section 14.11 we derive the Hamiltonian and discuss the spectrum. 
In section 14.21 we find that the wave function has an overlap with the ground state 
of the bilinear-biquadratic spin 1 chain close to unity in a region within the Haldane 
phase for = 6. In section 14.31 we compute the Renyi entropy and the two-point 
spin correlation function by Monte Carlo simulations. We find that the Renyi entropy 
grows logarithmically with the length of the considered subchain, while the two-point 
spin correlation function decays algebraically with the separation between the spins. 
The Hamiltonian for the chiral correlator of a product of spin 1/2 fields is given in 



Appendix B together with a discussion of the spectrum. Finally, concluding remarks 



are given in section [51 



2. Hamiltonians for the SU{2)k WZW model 

The WZW models, together with the minimal models of Belavin, Polyakov and 
Zamolodchikov, are the best known examples of rational conformal field theories 
(RCFT). In a RCFT the number of primary fields is finite, and they form a closed 
operator algebra characterized by the fusion rules. In this paper, we shall restrict 
ourselves to the WZW model based on the Kac- Moody (KM) algebra SU{2)k- There 
are k + 1 primary fields in this RCFT, which are denoted by (pj, where j = 0, |, . . . , | 
labels the total spin. Each field consists of 2j + 1 components corresponding to the 2j + 1 
possible values of the projection of the spin on the quantization axis. To refer explicitly 
to one of these components, we use the notation (j)jm, where m = —j, + ,j. The 
Virasoro central charge c, conformal weights hj and fusion rules are given by 

j=\h-h\ 

The simplest nontrivial WZW model corresponds to A; = 1 and has c = 1 and two 
primary fields 0o and 0i/2 with /iq = and /ii/2 = 1/4. This model can be easily 
constructed using a free massless boson whose central charge is c = 1. The next model 
in this series corresponds to /c = 2 and has c = 3/2 and three primary fields 00; 0i/2 
and 01 with = 0, hi/2 = 3/16 and hi = 1/2. This model can be constructed by 
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combining three Ising models (c = 3 x |), or alternatively it can be built using a free 
boson and an Ising model (c = 1 + ^). In this paper, we shall focus on the k — 1 and 2 
models, but the general formalism will be valid for any integer level k>l. 

2.1. Null vectors 

The rationality of the WZW models lies in the existence of null vectors in the 
representation spaces of KM algebras. The decoupling of these vectors provides the 
basic equations that allow for a complete solution of these models. Null vectors will be 
the key for the construction of Hamiltonians associated to WZW models, and for that 
reason we shall describe them in detail below. 

The SU{2)k WZW model is characterized by the chiral currents J" (2;) (a = 1, 2, 3), 
whose Laurent expansion J^{z) — Yl'^=-oo Jn^~^~^ contains the modes that satisfy 
the KM algebra 

[Jn^ >^m] = i XI ^"'"^ '^n+m + ^ab 5n+m,0, n, m G Z, (2) 

c 

where Sabc is the Levi-Civita symbol. Notice that the zero modes Jq form a closed SU (2) 
algebra. In the standard spin basis the algebra becomes 

[7-0 tO 1 _ ^„ r 

r jo j±i ^ ^ j± (3) 

where JO = J± = J^±iJ^. 

To construct the representation theory of this KM algebra, we first note that the 
primary fields introduced above effectively factorize into a product of a chiral field and 
an anti-chiral field. We shall only consider the former, which we denote by <i>j{z), where 
z is a complex parameter. In the context of conformal field theory, z represents a 
point in two-dimensional space-time, but later on, these coordinates will turn into free 
parameters in the Hamiltonians we construct. To each chiral primary field (pj{z), one 
associates a primary state satisfying the following conditions 

10,) = 0,(0)10), Jo"l0.)=^l0.), 4"I0.) = O, n>0, (4) 

where |0) is the vacuum of the theory (i.e., J">o|0) = 0) and arc the 2j + 1 
dimensional matrices of the spin j representation of SU{2) (the mth component of 

^10.) is E„.C„.'l0.w)). 

The representation spaces of the KM algebra are modules 1-Lj obtained by 
acting on the primary states with products of the negative mode operators, i.e. 
J"ni J"^2 ■ ■ ■ J-nM\^3) (^1) ^2, ■ ■ ■ , > 1)- These states are called descendants and 
can be organized in towers having a common value of the sum n — rii -\- . . . -\- tim that 
is called the Virasoro level, not to be confused with the level k of the KM algebra. 
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A null vector \x) is a descendant state which is also primary, i.e., satisfies the 
condition (recall @) 

J: \x) = 0, n>0. (5) 

Each module Hj contains infinitely many null vectors. Here, we shall consider the null 
vector [8] 

\Xj,j,) = {J-iT* = k + l- 2j, = k + 1- j, (6) 

which is the null vector occurring at the lowest Virasoro level. It is illustrative to verify 
([5]) for the case = . Using the commutator 

[Ji, [J^iY] = ^ {J^iY^' (fc + 1 - ^ - 2J°) , i> 1, (7) 
that follows from ([3]) by induction, one obtains, 

= {k + l- 2j) {r,f-'' (2j - 2J°)|0,,) = 0, 

where we have used that \(t>jj) is a primary state, Ji\(t>jj) = 0, and that its spin projection 
is j, JQ\<t>jj) = One can easily prove that 

■Jq I Xj,j, ) = 0, Jo I Xj,j. )=J*\ Xj,j. ) , (9) 

so that \Xjtjt) is the highest weight vector of a representation with spin j^,. In summary, 
the module Tij contains at level n^: a multiplet of 2j^: + 1 null vectors with spin j^:. 

Since a null vector is also a primary field, a null vector and its descendants comprise 
a module on their own, and this module can be shown to decouple from Tij [9]. The 
number of states dj^k{n) appearing at Virasoro level n after removing all the decoupled 
modules can be determined by considering the character 



oo 



n=0 

This character is given by [9] 

„{2j+l)2/[4(fe+2)] 

xAq, k) = ^- J2 [2^' + 1 + + 2)] q-^'^^'^'''~'^'^\ (11) 

n=—oo 

where 7]{q) is the Dedekin eta function 

oo 

V{q)=q'/''l[(^-qn- (12) 

n=l 

We collect the values of dj^k{n) for the lowest spin modules Tij and various values of n 
and k in table [TJ An inspection of the rows of this table shows that, for given n, the 
number of states remains constant beyond a certain value of k. This can be understood 
from the fact that states are removed whenever there is a null vector. For k < n* — l + 2j 
(recall ([6])), one or more null vectors are present for n < n*, whereas no null vectors are 
present for n < n* if k > n* — 1 + 2j . The number of states removed ai k = n* — 1 + 2j 
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Table 1. Number of states dj^k{n) for the modules Tij {j — 0,^,1). The value of 
dj,k{n*) appears with a star. The Virasoro level n labels the rows, and the KM level 
k labels the columns. 
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and n = n* is the number of null vectors in the multiplet |Xj,m*); i-e. 2j* + 1. This is 
expressed in the relation 

- C?i,n*-l+2j(^*) = 2j* + 1, Vfc > - 1 + 2j. (13) 

In the module Ho, for example, there are 4 states at n = 2 and k = 1, which turn into 
9 states for all > 1. The difference 9 — 4 = 5 is the number of null vectors that form 
the j^, = 2 spin multiplet at level n^, = 2 and k = 1. 

2.2. Clehsch-Gordan coefficients and projectors 

Equation ([6]) corresponds to the Clebsch-Gordan (CG) decomposition 

Vj, C \/®"* ® Vj (14) 

where Vj = (J = j^^l^j) is the vector space for the spin J irreducible 

representation (irrep). This decomposition consists in the symmetrized product of 
copies of the spin 1 irrep tensored with the spin j irrep to yield the highest possible spin 
j* = n^: + j. A basis of V^-, is given by (m^, = — j^,, — j^, + 1, . . . , j^,) 

\X,.m.)= Yl ^rX.,,.n^^\---^-rl0.™). (15) 
ai,...,an^ ,Tn 

The constants C*ai™*a„_^ can be constructed from the product of n^, CG coefficients 
and satisfy the condition 

^ai,...,an,,jm^ai,...,an:t,jm "mtin',^ \^^J 

ai,...,a„, ,m 

which ensures the orthonormality of the vectors Ixj.m,,)- We adopt the Condon-Shortley 
convention that all CG coefficients are real. We note that ( |T5l) could equally well be 
written in terms of the fields Xj*m,{^) 4'jmi.z), and we use this notation in the 
following. 

Let us next define the linear map 

Ct,: Vr'®V,^V,,, ^^^^ 
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which allows us to write fll6p as 

where l2j,+i is the identity matrix in Vj^. Multiplying C^^^j and (C:^*^^)"'" in reverse order 

Kn., = (C>,,.)'ct, (19) 
yields the map 

K„.j : V^f ® V^f ® Vj (20) 
that, using f lTSj) and f|T9|) . satisfies 

K^^,,. = K„,,,, = K„„„ (21) 

so that J- is a projector with components 

(IC \ai,...,an,,m _ \ ^ ^jtin^ /r,r)\ 

Using this equation, one can parameterize the null vectors f|T5|) with indices of the vector 
space Vj®"' ® Ij- as follows 

Xai,...,a„,,jm{^) = ^ ] a„, . jm Xj.m, (^) 

(23) 

bi,...,bn^,m' 

Since, by use of ( |T6|) . 

Xj,mt{^) = ^ ] ^ai,...,a„^ jm Xai ,...,a„,, ,jm(2^) ) (24) 

we see that Xa-i_,...,a„^,jm{z) fulfil the eigenvalue equation 

Xai,...,a„,,jm{z) = (■^"* j)6i,...,6„,*,m' Xfti jm' (^) • (25) 

,m' 

As a result of the above construction, K^^j must fulfil the following symmetry 
constraints 

/TV- ^...ai...aj...,m _ /-r^ y..aj...ai...,m _ /-p^ y..ai...aj...,m 

l-l^n*,jj...b;^...6,...,m' ~ l-^"*jV...6fc...bi...,m' ~ l"^". jV. ..?)(.. .bfc...,m' ' l^".* 

E (K„.,)C:;£:;;;:^ = O' E (Kn..)!:^:::::;::Z' ^m. = o, (27) 

''mm" l-^"*J/'6i,...,fe„,,m' ~ Z-^ K'^n* J J ...bi...,m' ^m'm" ~ \^°) 

ai,m" bi,m' 

and, indeed, ( I26ll28l) together with ( l2Ti) are sufficient to uniquely determine K„,j. The 
null vectors at levels k = 1 and 2 are collected in table [2], and one can check that the 
following ii'-tensors fulfil the conditions, 

{^2,o)llh2 = 2^^aibj_Sa2b2 + Saibi^aibi) - -Saiai^b^bi, (29) 



/ \ am 2 / ^ — -\ \ 

K^'Vbm' ^ 3 [ ~ ^ / abet mm' 1 



t^ = la", (30) 
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Table 2. Null vectors for fc = 1 and 2. The possible values of the indices are: 
a, ai, 02, as = 1, 2, 3, m = ±1/2, and m' — 0, ±1. 
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biM? = g i^aibiSa2b2Sa3bs + permutations of b's) - — , 
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b\b2m' 



3 1 

■ ('^aibi '^0262 ~^ '^ai b2 '^a2f'i ) 'Z^a\a2^b\b2 



■7 ^ ^ (^aibic (^0262 ^a2fe2c'^aifci ~l~ ^ai62C '^a2bi ~l~ ^a2&ic '^(11^2) ^mm' 1 



(31) 



(32) 



(33) 



where a" (a = 1,2,3) are the Pauh matrices. To simphfy the notation, we shall drop 
the m and m' indices in the following, considering that i^' is a matrix, and we shall let 
the summation over repeated spin component indices (e.g., 61, ... , &nj be implicit. 



2.3. Decoupling equations 

A chiral correlator is the vacuum expectation value of the product of a set of primary 
chiral fields: 

i){zi, ...^zn) = {(t)jAzi) ■ ■ ■ <PjN{zN))- (34) 

In our construction, ip{zi., . . . , z^q) is the wave function. Written more explicitly, the 
state of the spin chain is 

\ip{zi, . . . , zn)) = ?/'^j_...,^^(zi,...,zjv)|mi,...,mAr), (35) 

mi,...,m]v 

where the coefficients are given by 

'^mi...miv(^l5 • • • 5 ^n) = (</'jimi(^l) • • • '/'jjvmiv (^A^ ) ) • (36) 

We note that a chiral correlator is not necessarily one unique function. If the fusion 
rules ([T]) allow more than one possibility for how two fields can fuse (i.e., if there is more 
than one field on the right hand side of ([1])), there will be more different fusion channels, 
and each fusion channel corresponds to one wave function. All of these wave functions 
are ground states of the Hamiltonian that we shall construct. 

There are more approaches available to find explicit expressions for the chiral 
correlators flM|) . In some simple cases it is possible to solve the KZ equations |7|, 
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which are a system of partial differential equations in the coordinates Zi fulfilled 
by the chiral correlators. The KZ equations can be derived from the null vector 
X = — -^y^ J"]^ Jq)0, where L_i is one of the Virasoro generators. Another possibility 
is to use the Feigin- Fucks or Coulomb gas representation of the WZW model |9j . The aim 
of this section is to derive a set of algebraic equations satisfied by the chiral correlators 
that follow from the decoupling of null vectors in the modules Hj. These equations are 
algebraic, and not differential, as in the KZ case, and this fact will allow us to construct 
quantum spin Hamiltonians whose ground states are the chiral correlators of the WZW 
model. 

Before presenting this construction, let us notice that ( l3ll) is invariant under global 
SU{2) rotations generated by the total spin operator [8] 

N N 

T^'i^ = ■ ■■^jA^n)) = 0, T^^ (37) 

1=1 i=l 

where denotes the spin matrix acting at the site i, in the representation 
corresponding to the field 0j-(zj). We shall suppose for simplicity that ji = j, Vz, 
but the results can be easily generalized when the primary fields have different spins. 

Let us start by replacing the primary field (f)j{zi) in the correlator (IMIl by the null 
field (l23ll in the module Tij. One can show from the definition of a null field that a 
chiral correlator is zero if it contains a null field [9]. This is referred to as the decoupling 
of the null field. Therefore, 

{^,{z^) . . . x^'^'-'^-^'iz,) . . . (P,{zn)) = 0, (38) 

y„...,a„. ^ ^^^^^ ^^^^^ ^ ^ ^ ^ ^^^^ 

'"'/ bi,...,6„* 

for all Zi, . . . , zjv; where the superindex i in K^*^^^- implies that this matrix acts on the 
spin degrees of the operator 0j representing the ith spin in the chain. The modes J^^ 
act on the field located at the point Zj, but their action can be transferred to the primary 
fields located at the other positions by means of the Ward identity |8j 

(0,(zi)...(J^^)(z,)...0,(^;v))= ^^(0,(^i)---V^(^.)---0,(^7v)). (40) 

Here, and in all that follows, we use the convention that the sum over ii ^ 12 ^ ■ ■ ■ ^ 
inij^ in+ii ■ ■ ■ , in+m) mcaus that the indices outside (inside) the brackets should (should 
not) be summed over, and only terms for which all n + m indices are different contribute 
to the sum. Similarly, we write ii ^ 12 ^ ■ ■ ■ ^ in when all n indices are required to be 
different. Applying pUj) times to fl39|) yields 



I.e., 
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Vi;^;.---* {z,,..., zn) ^(^1, ...,zn) = 0, t = l,...,N, (42) 

where 

N 

Kr"* (^1, • • • , = E i^X'T '^■■■^A--- (43) 

V / b\...bn^ 

and Zik = Zi — Zk- These are the algebraic equations mentioned at the beginning of this 
subsection. 

Finally, we shall transform fj42|) into the following equation 

Ci^y^-' (^1, ...,zn) i;{z,, . . . , z^) = 0, t = l,...,N, (44) 

where 

N 

and Wik = (-2j + Zk)/{zi — Zk). To show this relation, let us apply the identity 

-^ = l+'-^ = l + w,k (46) 

to 

ai...a„. 2Zi \ - 2Zi K , / /.A 



T'^ - + E t^^nt-! I ■ ■ ■ ( - t?"* + E ^..„,tt: I ^ (47) 



where we have used (1281) in the final step. Next, we want to move the matrices to 
the right in this expression. This produces extra terms coming from the commutator 

[T^ ^ tl'] = ie,,,y^ E ""^^ ' (48) 

but the antisymmetry in the h indices causes these terms to cancel out when multiplied 
by the tensor K whose indices are totally symmetric. Once the T^^ matrices are brought 
to the right, they annihilate the correlator if), which proves that 

(2^,)"* v';y-' = dil]--"-* (49) 

and thus (jHj) is equivalent to fH2l) . 
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2.4- Spin chain Hamiltonians 

We can use the operators in ( l45l) to define the matrices 

= {ci:r"* (^1' • • • ' ^^)) ' c!^-"* (^1' • • • ' ^^)' (50) 

which are Hermitian, positive semidefinite and rotationally invariant 

= <,>0, [T^<,]=0. (51) 

They annihilate the correlator ip, 

H^:1^^{zu...,zn) = 0, (52) 

and their sum can be taken as the Hamiltonian of a spin chain with spin j at each site, 
i.e., 

N 

H^„ = Y,Htr (53) 

i=l 

It follows immediately that is the ground state of this Hamiltonian. If the fusion rules 
give rise to several chiral correlators, then i^n* ,j will have several ground states in one-to- 
one correspondence with the chiral correlators. Finally, we note that Hamiltonians that 
are translationally invariant by one site can be achieved by choosing the coordinates Zk 
to be equally spaced on the unit circle, i.e., 

Zk = e^'^^^/^ {k = 1, . . . , N) =^ Hn,,j : translationally invariant. (54) 

We shall refer to this choice as the uniform case. 



2.5. Algebraic equations for spin correlators 



The decoupling equations (jS]) also imply certain relations between the different n-point 
spin correlators. Let us write 



c 



«,ai...a„, 

n*,j 



Arr"*(^i'---'^iv)+Qr: 



as the sum of an anti-Hermitian and a Hermitian operator. 



ii,ai...a„^ 



and define 



T 



'bl...bn 



t'^ 



ii,ai...a„^ 



j,ai...a„. 



t 



-ll...ln ll 

The decouphng equations (l44l) imply 
and adding and subtracting these equations, we find 



pi,a\...an 



T 



•bi...bn 



+ 



rpb\...bn 



j,ai...a„. 



.j,ai...a„. 



_|_ J rpbl...bn 0i,CLl---Cinf 



0, 
0. 



zn) (55) 
(56) 
(57) 

(58) 

(59) 
(60) 



actual expres 

of n, (l59l) and (1601) provide a set of linear equations between different spin correlators 



Inserting the actual expressions for All'^j""'"* and Q.^"^"""* and choosing different values 
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3. Results for the SU{2) WZW model at level k = 1 

We now turn to a more detailed investigation of the case k = 1. First, we provide explicit 
expressions for the wave function and the Hamiltonian for N spin 1/2 fields. We then 
show that the model is a generalization of the HS model and discuss the degeneracies of 
the spectrum. As mentioned in the Introduction, the overlap with other spin models, the 
Renyi entropy and the two-point spin correlator have already been computed numerically 
in [6], and consequently we here concentrate on deriving analytical and semi-analytical 
expressions for the spin correlators, which is done in the last part of this section. 



3.1. Wave function and Hamiltonian 

At level k = 1 there are two primary fields 0o and 0i/2 whose fusion rules are 

01/2 ® 01/2 = 00, 01/2 ® 00 = 01/2, 00® 00 = 00- (61) 

These rules tell us the number of chiral correlators that arise from a product of primary 
fields on a given geometry, which in our case is the sphere. For k = 1, the nontrivial 
correlators are obtained by multiplying N fields 0i/2- Equation fl6T|) implies that the 
fusion of these fields produces the field 0o if N is even and 0i/2 if is odd. Since the 
vacuum expectation value of 0i/2 is zero, and that of 0o is not, A^ must be even to get 
a nonvanishing correlator 

^si,...,sjv(2;i, ...,Zn) = (01/2,si(2i) • ■ ■(Pi/2,sn{zn))- (62) 

Here, we have made a slight change of notation compared to the previous section, such 
that Si/2 = ±1/2 is the third component of the i^^ spin. In addition, '0si,.--,siv(^i, • • • , ^n) 
can only be nonzero for YliLi Si = because -0 is a singlet. Note that the correlator 
( l62l) is unique since the fusion decomposition only generates one field at a time. 

The explicit form of ( l62l) can easily be obtained from the free field representation 
of the SU{2)i WZW model in terms of a free boson (p{z) whose two-point correlator is 
{ip{zi)ip{z2)) = — ln{zi — Z2). In the Cartan-Weyl representation, the current algebra 
operators J"'{z) are given by 

J\z) = ^d,^, J±(^)=e±^v^^(^), (63) 
and the primary fields 0i/2,s, (si = ±1) by 



where : . . . : denotes normal ordering. The two-point correlators of these fields are 

{J\z,)j\z,)) = {J^z,)J-{z,)) = (65) 

^. ^isM^i)/V2 .. ^isM-2)/V2 ^ ^^^^^^^^ ^.1^2/2^ ^gg^ 
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where 2:12 = Zi — Z2 as above. Plugging into ( 1^ . one gets [H] 

N 



13 



(67) 



i<j 

N : even, 



AT 

E 



0. 



The Marshall sign factor pi ensures that ip is a. singlet state, as required by (I3j 

The Hamiltonian H^ i, whose ground state is given by f l67|) . can be constructed 
using the results of the previous section. For \zi\ = 1, Wi, we find 



N 



(68) 

+ (69) 



AT 



4-Hi/2 = ^1 1 



AT 



TV 



E4H-E 



^ Wfci t," t"j 



■ (70) 



3.2. Relation to the Haldane-Shastry model 



In 1988, Haldane and Shastry considered the state obtained by the Gutzwiller projection 
of the one dimensional Fermi state at half filling [21 S] (see [TU] for a review) 

I^G) = n(l-^^.t^a)|FS), (71) 

i 

where n^^s = cj^Q^s (s =t, i-,i = 1, • • • , -?V) are the occupation operators of an electron in 
a periodic chain with N sites, and |FS) = n|/c|<A:^ the Fermi state at half filling. 



The state (I7T1) has a unique spin per site because the charge degrees of freedom are 
frozen by the Gutzwiller projection. These authors then wrote fl7T|) in terms of hard- 
core boson variables. Mapping the spin up (down) into the empty (occupied) states of 
a hard-core boson, they found 

m= Yl Hnu...,n^/2)bl...bi^^j0), (72) 

ni<...<N/2 

where fej^ is the hard-core boson creation operator at the n^^ site and 

2 



^{ni, . . .,nN/2) oc e 



N/2 
ni>nj 



sm 



TT{ ni - Uj) 

N 



(73) 



Using this equation, Haldane and Shastry showed that \iPg) is the ground state of the 
spin Hamiltonian 

Zi Z'j 



HS 



(74) 
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where Vij is the permutation operator of the spins at the positions i and j. The 
permutation can be expressed in terms of spin operators as Vij = 2t°;t'^ + 1/2, so fTMl) 
takes the form 

1~ 



H 



HS 



ZiZj 



2t'^t] 



(75) 



In [H], it was shown that the wave function fl75]) corresponds to the state fl^7|) for the 
choice Zn = q'^™/^ , Indeed, let us denote by = 0, 1 the occupation number of a 
hard-core boson at the site i. The mapping between the spin variable Sj = ±1 and is 
given by Sj = 1 — 2qi and one can easily show that fl73l) becomes 



ip{si, . . . , sn) oc e'2 E«:odd*> 



n>m 



iT(n — m) 



sm 



N 



J]s, = 0,(76) 



where the constraint J^iLi = guarantees that the total third component of the spin 
vanishes. This expression is proportional to fj67|) . 

From these considerations one must expect the HS Hamiltonian ([71]) to be 
intimately related to the Hamiltonian (!70|) . To show this relation in detail, we need 
some simple identities. For generic values of Zi, one can show that 

.2 



1 + 4 



ZiZj 



and 



WkiWkj + WijWik + WjkWji 

The latter identity implies 



N 



N 



Wki Wkj = 2Wij + Wij{ci - Cj) + N -2, Q = ^ 



Wki, 



which plugged into (!70|) yields 

^1/2 = 



1^£^ iV_ (iV+1) 
9 r2. ^ 4 



2 "7" 



6 



W^^j ( C'i Cj ) 



where is the total spin operator defined in ( |371) . In the uniform case, we have 



Vi, 



N{N^ - 1) 
12 ' 



so 



^1/2 — ^ 



HS 



-En 



En 



N^ + 2N 



(77) 



(7^ 



(79) 



(80) 



(81) 



2) 



6 ' - 12 ' 

where £^0 is the ground state energy of "Hhs- The positivity of 7^1/2 implies a lower 
bound for the energies of the eigenstates of "Hhs with total spin S, i.e.. 



Es-Eo>-{N + l)S{S + l), 




(83) 
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which is indeed satisfied by the lowest eigenenergy Es in the total spin sector S < N/2, 



Eq — En 



-S(2 + 3NS -2S^). 
6 



The latter formula can be derived from the general expression for the eigenvalues of 
( FM|) . The eigenvalues were found in [11] and can be expressed in terms of a set {rrij} of 
nonconsecutive integers rrij G {1, 2, . . . , — 1}, nij+i > rrij + 1, called rapidities [I2] 

Ei{m,}) = J2rnj{m,-N). (85) 

j 

Thus, for example, the ground state corresponds to the set {rnj} = (1, 3, 5, . . . , — 1) 
(iV even), and the lowest excited state of total spin 1 to {rrij} = (2, 4, 6, . . . , — 2). 

In the nonuniform case, where Cj 7^ for some i, we get a generalization "Hhs of the 
HS model. Specifically, 

{N + 1 



6 



' rpar-pa 



where 



HS 



E 



Cj) 



and the ground state energy is 



+ 



24 



N 



6) 



(87) 



(88) 



Note that the coupling strength in flHTI) is more complicated than just the inverse of 
the square of the distance between the spins measured along the chord. In p], a 
different approach based on the KZ equation also led to the Hamiltonian "Hhs; cind 
it was demonstrated that ip is an eigenstate of Tins with eigenvalue Eq. We provide the 



details of this derivation in Appendix A 



The spectrum of the uniform HS model exhibits an enhanced symmetry given by 
a Yangian Hopf algebra [IB] generated by the total spin operator and the operator 

m 



A" = ^^WijeabctWj- 



The operators T" and A" all commute with the Hamiltonian fTMj) but they do not 
commute among themselves, and they are the key to compute the complete spectrum 
mentioned above. Observe that A can be constructed from the operator ( 168|1 . 

lE^n- (90) 



It follows that ijj is an eigenstate of ( 190|) (with eigenvalue zero) in both the uniform 
and the nonuniform cases. In the nonuniform case, however, the Yangian symmetry is 
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broken, meaning that (1901) does not commute with (18 7p . In the last part of Appendix A 
we show that is also an eigenstate (with eigenvalue zero) of the operator 



( I ^ 17 



which is a generalization of the Inozemtsev invariant to the nonuniform case. In 
the uniform case, this operator commutes with the Hamiltonian (17il) [S], but does 
not commute with "Hhs in general. 

To illustrate the breaking of symmetries, we compute the spectrum for the three 
different choices of Zi considered in [6]: uniform, dimer and random. The dimer model 
is defined by the following choice of parameters 

e2.ib+(-ip5]/iv^ 0<5<-, J = 1,2..., AT, (92) 

where 5 characterizes the dimerization of the ground state. The value of Cj defined in 
(I79D is 

Cj. =i^ (-1)^+1 tan(7r(5), (93) 

which only vanishes in the uniform case 5 = 0. The random model is defined by the 
choice 

5 . 5 



^2.iU+^,)/N ^ -2<0. <2' J = l,2...,N, (94) 



where 0j is a random variable uniformly distributed in the interval [— |, |]. The 
parameter 6 determines the degree of disorder in the ensemble. The choice ( l94l) gives 
the random HS version of the spin 1/2 antiferromagnetic Heisenberg model (AFH) with 
random exchange couplings. The latter model belongs to the universality class known 
as the random singlet phase. The Renyi entropies of this phase differ from that of the 
uniform AFH model given by the SU {2)i model [151 116]- 

The spectrum of the generalized HS Hamiltonian ( 187|) for these three models is given 
in table |3]for A^ = 6. The Hamiltonian "Hhs; the dot product T'^T'* of the total spin 
operator with itself and the z-component T° of the total spin operator all commute with 
each other. It is hence possible to diagonalize these operators simultaneously and label 
the resulting states by their energy E, spin S and z-component of the total spin Sz- We 
have omitted from the table, because also commute with the Hamiltonian, and 
it follows that the 25* + 1 states within each spin multiplet are degenerate. Instead we 
write only the number of states in each multiplet. In the uniform case, the momentum 
operator also commutes with T" and the Hamiltonian due to invariance under translation 
by one site. The momentum operator can be defined as pu = — iln(T), where 

T = Vn,N-iPn-1,N~2 ■ ■ ■ "^2,1 (95) 

is the translation operator and Vij is the operator that permutes spins i and j as defined 
just after ( |74l) . The momentum eigenvalues in units of 7r/3 are given in the table. The 
fact that the chain is invariant under the flip operation 

'Pn,i'Pn-i,2 ■ ■ ■ 'Pn/2+1,N/2 (96) 
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Table 3. Spectrum of the generalized HS Hamiltonian (|87l) for = 6 sites and three 
different choices of zj: Left: uniform ((8T|). Middle: dimer with 6 ^ 0.1 ((92|). Right: 
random for one particular realization of and (5 = 1 (|94p. E is the energy, S* is the 
spin, p„ and pd are momenta (see text), and # = 25"+ 1 is the number of (degenerate) 
states within each of the spin multiplets. 



E 


s 


3Pu 

TT 


# 

7/ 


E 


S 


3pd 
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IT 
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3 


-16.3598 
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-2 
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-2 
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-2 
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1 
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3 
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1 
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-11.4068 
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1 


3 
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-8.9570 
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-8 
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-2 
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-8.8204 
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-2 
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-7.9934 
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-2 
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-7.7381 
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-7.5353 
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-5 
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-1 
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-2 
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-5.0000 
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-4.8865 
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5 





3 
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-0.0000 


3 





7 


-0.0000 


3 


7 



explains why states with momentum eigenvalue ±7r/3 or ±27r/3 always appear in 
degenerate pairs. In the dimer case, Pu does not commute with the Hamiltonian, but 
since the chain is invariant under translation by two sites, we can instead add the 
momentum operator p^ = — iln(7~^) to the set of mutually commuting operators. Also 
in this case, states with pd eigenvalue ±27r/3 always appear in degenerate pairs due to 
invariance under the flip operation fl96l) . For the random case, there is no translational 
symmetry. Looking at the degeneracies between different spin multiplets, we observe 
that the number of degeneracies decreases in going from the uniform to the dimer case 
and in going from the dimer to the random case. In the dimer case, the Yangian 
symmetry is broken and only degeneracies between pairs with opposite pd eigenvalues 
occur. In the random case, there are even less symmetries, and all spin multiplets have 
different energies. 

This subsection shows that the generalization of the uniform HS model embraces 
several interesting models within a common framework. We shall show below that 
more information about this model can be obtained exploiting the algebraic decoupling 
equations. 
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3.3. Spin correlators 



In 1987, Gebhard and VoUhardt computed the two-point spin correlator in the 
Gutzwiller state ( 17T|) for an infinite chain, obtaining the exact formula [IT] 



Si(7™) 
Ann 



'"n 2 1^,3 ^ s 



where Si{z) is the sine integral function 

sin(t) 



dt 



t 



(97) 



(9^ 



The appearance of the term 6ab refiects the SU{2) symmetry of the Gutzwiller state. It 
may be worthwhile to notice that this result triggered the Haldane-Shastry works [HIH]. 
Using the asymptotic expression of the sine integral function 

_l)n+l 



one obtains 



Si(7rn) 



nn 



n > 1, 



1 



6n 



> 1. 



(99) 



(100) 



One can also derive fllOOp from the underlying CFT that describes the long distance 
behaviour of the HS model that is given by the SU{2)i WZW model. The spin 
operators can be expressed in terms of the current operators J"'{z) (!63l) and their 
antiholomorphic version J'^(^) as [T8| 119] 



^ J'^iz^) + J%z^) + (-1)^ 



IT 

— COS 

2 

TT 



J^{Zn) + .r{Zn) + (-l)%/-exp ±i 



V2 



(101) 



where z„ = 2TT'm and ip{z,z) = ip{z) + ip{z). It is interesting to compare the spin 
correlator of the HS model with that of the AFH model. The latter correlator behaves 
as [SniEl] 

-l)'^(ln(n))i/2 1 



c- 



n 



(102) 



where C is a nonuniversal constant. The logarithmic correction to the leading 1/n 
term is due to the marginal irrelevant operator J""]"", which added to the action of 
the WZW model gives the long distance behaviour of the AFH model [H]. In this 
respect, the HS model is at the fixed point of the renormalization group, so that the 
continuum Hamiltonian does not contain the marginal operator J" J", which explains 
why the logarithmic corrections in the correlator fllOOp are absent. 
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3.3.1. Algebraic equations for spin correlators We shall next derive f l97p and several 
other expressions for spin correlators from f l59|) . From f l68|) . we find 

where it is assumed that \zi\ = 1, such that Wij is imaginary. Note that A^'^ and Q''*^ 

-•^'2 '2 

are one-body and two-body operators, respectively. For the spin 1/2 product 
of spin operators acting on the same spin can always be written as a sum of terms 
containing at most one spin operator each, so it is sufficient to consider spin correlators 
for which all the spin operators act on different spins. We hence assume that all the 
indices ii . . . in T°-^--f" are different and also different from the index i. T"-^-f" jg then 
an n-body operator. If we choose n odd in fl5^ . we get 3"'+^A^(A^ — 1) ■ ■ ■ (A^ — n) linear 
equations involving only (n + l)-point and (n — l)-point spin correlators. Starting from 
n = 1, this allows us to successively compute spin correlators with higher and higher 
numbers of operators as we shall see below. 
Finally, we mention the important result 

which follows from (159!) with n = 1 and bi = ai. (11041) is a system of linear equations 
for the two-point correlators with i ^ j that can be solved numerically for 

generic values of the parameters zj, and analytically in special cases. In the following 
paragraphs, we shall investigate some particular solutions. 

3.3.2. Two-point spin correlators for the infinite chain Let us choose 

^. = e2--./^, x,-GM, j = -^ + l,...,^. (105) 
Then, in the limit 

l^il < iV, V«, =^ Wij = -icot ^ ^ , (106) 



•^3 



(I104p becomes 



where the equality is reached in the limit — )■ oo. In the uniform case, i.e. 
x„ = n (n G Z), translational invariance implies that Cij only depends on the difference 
\i — J I, hence it is sufficient to consider the case z = and j = n in (I107p . 

C(n)+ ^ ^C(n-fc) = 0, n^O, C{n) = Co,n = 0^,0. (108) 

where we have used that C(0) = {f^t^) = 3/4. To solve these equations, we introduce 
the Fourier transform of C{n), 

G{q) = e'-'^Cin), C{n) = ^e-'"^G'(g). (109) 
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Using the identities 

J2 ne^"^ = -27ri-5(g), ^ -e"^ = ivr (sign(g) - ^) , (110) 

n=— oo n=— oo 

fllOSp turns into the differential equation 
/ . q\ dG{q) 3 

The SU{2) invariance of the ground state implies that G'(O) = {T"-tQ) = 0, hence 
integrating (11 111) yields 

G(g) = -^lnfl-My (112) 



4 V ^ . 

Finally, C{n) is given by the inverse Fourier transform, 

C(n) = r ^e-"« In fl - (113) 



4 27r V ^ . 

— [ do cos(nq) In (l — — 

3 r 1 • . N 1 3 , ^ r" , sin(t) 

/ dg sm(ng) — = cos(7m) / at . 

Jo 1 - g/TT 4?™ Jq t 



Since (t^tg) = Sabit^t^o)/^ due to the SU{2) symmetry, ([97D follows. 

3.3.3. Two-point spin correlators for finite chains Let us next show that (I104p can be 
solved analytically for a uniform spin chain with a finite number of sites A^, and not 
only in the large N limit, as done above. We write 

zj = z^ (j = l,...,iV), ^^e^'^^/^ (114) 

so that (I104p becomes 

^ k=l ^ 

If n = N/2, the first term of this equation vanishes, which leads us to split this equation 

as 

1 -n 1 _|_ ^-k j\T 

^H + E i + ,-niI,-. ^("-fe) = Q- -^Y'^' (116) 
^-1 1 I -k 

E I3^^W2 - ^) = 0. (117) 

k=l 

To solve these equations, one introduces the finite Fourier transform (recall (11091) ) 

N N 

G, = Y.z-^C{n), C{n) = -J^ (l^^) 

n=l g=l 
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where Gj\f = J2n ^(^) ~ 0- Now we multiply f lll6p by 2;"'^ and sum over n 7^ N/2, N to 
obtain 

G,-l-(-irC(N/2) + ^^G, Y: J- ,vi±i_^ 0.(119) 

q'=l n=l{^N/2) k=l 

while (I117P becomes 



5:g,,(-i)'X-"'1^=o- (120) 

(j'=l A:=l 

The two sums appearing in (11191) are the discrete version of the integrals flllOp . and 
they are given by 

^-1 , , _fe 



^(9)-ivS"''T^ = ' + ^-^-77' 1^^^^' (121) 

k=l 

B{Q)^l^ E -"^1^ = (-1)^ Gig-(^) - I) ' k\<N-h (122) 

n=l{^N/2) ^ ^ 

in terms of which (I119p and (I120p become 

- - + (-l)^a + i-iyNj2G,'A,i-ir' sign(g - g') = 0, (123) 



4 

g'=l 



AT 



J2G,'AA-iy' = 0, (124) 

g'=l 

where a is a constant given by 

N 

a = 2j2GqAgq{-iy -C{N/2). (125) 

9=1 

It is convenient to define the quantity 

gq = {-iyGq{q = l,...,N), (126) 
which must satisfy 

3 



gg - -i-iy + a + Nj2 sign(g - q') = 0, (127) 

q' = l 

N 

Y^gq,Ag, = 0. (128) 

Taking the difference for two consecutive values q and g + 1 in (11270 yields 

gg+i -gg + N{gqAq + gg+iA^+i) + ^(-1)« = 0, 1 < g < iV - 1, (129) 
which, when translated into the Gq quantities and using fll2ip for Aq, becomes 

(l - ^) iGq-,1 -Gq) = ^, 1 < g < iV - 1, (130) 
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Figure 1. Two-point spin correlation function C(n) = {f^tg) (n = 1, . . . , — 1) for 
the uniform HS wave function for a chain with A'^ = 20 sites. The results correspond to 
(i) the exact formula (jl32l) . (ii) the antiferromagnetic Heisenberg model with nearest 
neighbour couplings, and (iii) the thermodynamic limit (j97p . 



and hence 



N-l 



2^ N-l 

q'=q 



1 N 



G 



N 



0. 



(131) 



One can check that (11251) is also satisfied and that the parameter a is equal to 3/4. 
Plugging this expression into (11181) . one finally finds the correlator 



{n) = ^8a, C{n) 



6, 



ab 



-1)' 



4A^sin {nn/N) 



N/2 

E 

m=l 



sin 



m- 2 



(132) 



In the limit — > oo, one recovers fl97j) . Moreover, one can take the limit ^ 1 and 
n ^ 1 with the ratio n/N kept constant, obtaining 

7r(-l)" 1 



(133) 



]N sin (nn/N) AN^ sin^ {Tin/N) _ 

Equation (I133P can be obtained from the bosonization formulas (llOip with the 
operators J"'{z), J"'{z) and cos((y9(z, z)/\/2) being defined on a cylinder of circumference 
27riV, i.e., z = z + 2'KiN . The two-point correlators needed to reproduce (I133P can be 
deduced using the conformal transformation from this cylinder to the complex plane, 
z ^ w = e^/^ , which yields 



(J°(27rini) J°(27rin2))cyi ■ 
if{2'K'mi, — 27rir?,i) 



8A^^ sin^ 



cos 



V2 



cos 



(y9(27rin2, —2Tim2] 

71 



cyl 



4A^ 



sm 



-K^rii - n2) 
N 



-1 
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Figure 2. Absolute value of the two-point spin correlation function (i.e., (— l)"(t$5,to)) 
for the dimer model with (starting from above) ^ = 0, 0.01, 0.02, . . ., 0.10 and N = 4000 
(only the region 1 < n < N/2 is shown). The dotted lines are all proportional to 

[sin(7rn/iV)]-2 



As a summary of these results, we present a comparison between the correlators 
computed with the formulas f l97|) and fll32p for the HS model and the correlation function 
for the AFH model in figure [H 

3.3.4- Two-point spin correlators in dimerized chains In the dimer case ( l92l) . 
translational invariance by one site is broken, but there is still translational invariance 
for translation by an even number of sites. From the symmetries of the wave function, 
it follows that (t^t^) = (tm-nPo) ^ odd and (t^t^) = i^n-m^o) even, so 

it is sufficient to compute (t^to) — 1,2, ...,A^. Figure |2] provides numerical 

results for this correlator for a chain with = 4000 sites obtained from f ll04p . For 
5 = 0, we observe that the absolute value of the correlator decays as in the region 
I n N as it should according to f llOOp . For 6 = 0.1, we find that the absolute value 
of the correlator decays as [sin(7r?T,/A^)]~^, when the spins are sufficiently far apart. 
As 6 decreases towards zero, n/N needs to be closer and closer to one half for this 
approximation to hold. It is expected that the decay is algebraic for 1 <^ n <^ A^ 
rather than exponential even if the model is not critical because there are long-range 
interactions in the Hamiltonian [22| 123] . 

3.3.5. Four-point spin correlators We shall now use ( l59ll to compute the four-point 
spin correlator {t1^ f^^ 1°^^ t"*). As explained above, we can assume the indices zi, . . . , Z4 
to be different and consider the quantity 




(134) 
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Choosing n = 3 and T3 = t'^H'^H'^i^ (ii,i2,«3 all different) in yields 

0= ^ 2 Wj^fc Cj"^-"^"^"'' + -Wj4i^5aia4<^j"^^^^ (135) 

fc(7^j4) 

(s /^ba2azb /~ia4,a2as,ai\ _i_ (x /^aiba:ib f^aia4a:ia2\ 

^aia2bb _ ^01020403 N 
"r"'«4«3 I,"a3a4'-"iij2«3«4 n«2*3«4 / ' 

where Cf^'j"' = {ti^f^^) = 6aia2Ci% = ^aia2C{\ii - i2|)/3 is the two-point spin correlator 
found above. 

The expression can be simplified noting that, by rotational invariance of the ground 
state, one has the decomposition 

^^01020304 _ r r ^(1) I c r ^(2) , r r ^(3) ( T QfW 

^ili2hiA ~ "aia4"a2a3'-^iij2«3j4 ^ "a2a4"aia3 '-"11121344 ^ "a3a4"aia2 '-"iii2i3j4 ' 1"^"^^^ 

where C\'^]^^,^^_^ (a = 1, 2, 3) are invariant tensors that can be obtained from Cf^^j^-^j^"" as 

/^(l) ( A/^abba /^abab /^aabb \ 

nj2«3«4 ~ 20 V n«2i3«4 «ii2i3«4 hhhu) ' 

/^(2) 1 / /-labba _i_ A/~iabab /-~iaabb \ (^'\'7\ 

iii2i'iii 30 *i*2«3«4 ~^ iMisii iMiaU) ' \^'^') 

^(3) / /-^abba /^abab , A^aabb \ 

ili2i3ii ~ 20 V «l«2«3«4 ~ '-"iii2«3«4 ^iii2hU) ■ 

Permuting the spin matrices in (11341) . one can easily relate these tensors, 

'-'n«2*3«4 — ili2*3*4 ni2*4*3 «li4i3j2' V"'^'-'"/ 

r' . . . = r* . . . . = r* . . . . = r' . . . . 

^«l«2«3«4 '-'12«1«423 ^13«4«1«2 "-"«4«3«2«1 • 

Plugging f ll36p into f ll35p . one finds an expression, where the terms proportional to 
^0104^0203, etc, must cancel independently. This leads to 

0= J2 '^^^.kC^^2^sk + ^^^.nCni2^^U + (^un-W^,^^C^^^^^^^ (139) 

,(3) , «^i^{0) 

2 

0= Yl 2«^M;^C^^i.3/c + 2^M.2CS.3H + K4.2-«^Mn)cS^^^^ (140) 

^(y,. . _ y,. . \r^^) + ^^4^2 ^(0) 

I "'14*2 "'«4«3/'-"jij2j3i4 ^ 2 «1«3' 

0=^2 Wi,k Cl^^^^ + 2wi,i^Cl%i^i^ + {wi,i^ - Wi,i,)C^lli^j^ (141) 

If,,,. . ^r^^) + ^^Mo) 

These equations can be solved numerically for the four-point correlators in terms of the 
two-point correlators which are known from the solutions of (11041) . We can also exploit 
the relations (I138P to rewrite the above system in the compact form 



T^V"'l4il ^tAt2)'^ni2i-j,i4, ^ 9 ^j2«3' 



«4n ^(0) 



+ (w^i4il Wij4i2)C*iij4i3i2 "I" Q C*, 



«2«3 ■ 
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Figure 3. The four-point correlation function (iii5o^i25^n) in the uniform model 
computed from the exact formula (I142p and the approximation (11441) . which is valid 
when all the spins included in the four-point correlator are separated by many sites. 



For the uniform model, we can compare these formulas to the CFT prediction 
for the correlator of M spins. We shall only consider the most relevant term in the 
expression ( llOip for namely 



+0 _ 

''n — 



71 ( ip{Zn, Zn) 
— COS 1= — 

2 V V2 



(143) 



where Zn = 2nin {n = 1,...,N). Employing the aforementioned map between the 
cylinder and the plane one can derive 



\j=i 



V2 



(144) 



cyl 
M 



sm 



71 {Ui - Uj) 



Si,...,SAf=±l Kj '- 

A comparison of this formula for M = 4 and the results obtained from the numerical 
solution of f ll42p in the uniform case is given in figure [31 



4. Results for the SU{2) WZW model at level k = 2 

The SU{2)2 WZW model has three primary fields (pj associated to the spins j = 0, 1/2 
and 1. The nontrivial parts of the fusion rules are given by 

01/2 ® 01/2 = 00 © 01, 01 ® 01 = 00, 01/2 ® 01 = 01/2- (145) 

This theory is equivalent to the product of the critical Ising model and a free massless 
boson, so that their Virasoro central charges add uptoc=l/2-|-l = 3/2. The current 
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operators can be written in terms of the Majorana field x{^) of the Ising model and 
vertex operators of the chiral boson as 

J°(z) = iOM^), J^{z) = V2 x{z) e^''P^'\ (146) 

In the following, we investigate the model obtained from a product of spin 1 fields, 
and we note that the primary fields with spin j = 1 and conformal weight hi = 1/2 can 
be constructed as 

/:e±^^(^»): = ±1 _/ 1 ^ : even 

[ Pi,iX{^i) -^i = -1 2 : odd 

For brevity, we restrict the numerical analyses to the uniform case in this section. The 
Hamiltonian for the model obtained from a product of spin 1/2 fields (/c = 2, j = 1/2) 



is given in Appendix B 



4.I. Wave function and Hamiltonian for j = 1 

The above representation of the primary fields enables us to compute the chiral correlator 
of A^ spin 1 fields with spin projections Si = ±1, and parameters [5l [9]: 

^.,...,.^(^1, . . . , ^iv) = Pi(: e^^^^(^^) : . . . : e^-^^^-) :)(x(^i)'-^' • • • xM'-'h 
^ / 1 \ 

Hiz, - z,r'^ Pfo ^— , (148) 



Pi 



Zi Zj 



N 



pi = (-l)^-°dd{''» A^ : even, = 0. 

i=l 

The Pfaffian Pfo(. • •) comes from the correlator of the xi^) fields and is hence restricted 
to the positions where = 0. The sign factor pi ensures that ip is a singlet. The wave 
function vanishes if A^ is odd because an odd number of (pi fields fuse to a single 0i 
field. We note that a different wave function containing some of the same elements was 
investigated in [21] , where a three-body Hamiltonian, whose ground state is close to the 
wave function, was found numerically. 

To find the explicit expression for the Hamiltonian, we rewrite fl5^ into 

^Ai - ^i'J.s^ - + «'*?). (i-jfc = -iEafc, (149) 



(a1;i 



b 3 12 ' 12 



and consequently 

Jl r9 ^ . 1 . . .1 

(150) 



^ -2 , 5 . , ^ 1 , , , , 



Ci^i = ^ Wik 

When we plug this into (!50|) . assume = 1, Vi, and exploit the relations 



\t'i-,ti] — ''^^abct'i-, — 2, Eabctiti — iti, (151) 
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Table 4. Spectrum of the spin 1 Hamiltonian ()153p for Zj — ex-p{2TTij / N) and = 6 
(only the lowest 106 states are given). E is the energy, S is the spin, pu is the 
momentum, and =ff = 2S + 1. 



E 


S 


3p.. 


# 


table continued. . 




table continued. . 
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3 
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3 
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2 


-2 
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3 
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7 


25.5448 
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-1 
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we find 



N ^ N ^ N ^ N 



Hfl = -\Y.<-\ E ^^^^^^n-\Y.^m + \Y.<^^ (152) 

1 ^ 

+ g 5Z W,jW.,kttt''fit\. 

The Hamiltonian for the spin 1 case is thus 

N N / N \ N 

Hi,i = E 4 - ^ E h^- + 2 E + ^ E ^^'^t^)" (153) 

1 ^ 

+- Yl w^.w.utlt'^t^tl. 

The spectrum of (11531) for the uniform case and = 6 is given in table HI Since Hi i 
commutes with T", we again have that the 25 + 1 states within each spin multiplet 
are degenerate. Furthermore, we can also classify the states by the eigenvalue of the 
momentum operator pu = — ilii(T), where T is the translation operator defined in ( l95i) 
with the permutation operator given in the spin 1 case by [25] 



'Pid = (ttf^y + ttt] - 1. (154) 

We observe that all degeneracies between multiplets can be explained from the symmetry 
under translation and under a flip of the spin chain. This suggests that there is no 
Yangian symmetry in this model. We have tried to restore the Yangian symmetry by 
creating accidental degeneracies in the spectrum by giving the terms in the Hamiltonian 
(|53|) different weights, but we did not find a solution. 

4-2. Comparison to the bilinear-biquadratic spin 1 chain 

In [6], it was suggested that the dimer version of the generalized HS model for suitably 
chosen 5 could be close to the spontaneously dimerized phase of the J1-J2 Heisenberg 
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Figure 4. Overlap KV'IV'^)! between the ground state ip of (|153p and the ground 
state ipif, of (|155p as a function of for iV = 6 sites (when the ground state of 
is degenerate, we choose the state in the ground state subspace, which maximizes 
1(^1^^) I). The dotted vertical lines indicate the transitions between different phases 
of the bilinear-biquadratic spin 1 chain in the thermodynamic limit. 



model because the Hamiltonians of both models involve terms of the form and 
because tuning 6 can change the pattern of coupling strengths. It turned 
out that there indeed is a significant overlap between the wave functions of the two 
models, and one may ask whether the k = 2 and j = 1 model resembles another model. 
In the uniform natural choice for comparison is the bilinear-biquadratic spin 1 

Hamiltonian [26l EH [28] 



with periodic boundary conditions = t^. We compare the state ip in (11481) to the 
ground state ip<j> of ( 1155^ for = 6 in figure HI The phases of the bihnear-biquadratic 
spin 1 model for N ^ oo are also indicated in the figure. In the Haldane phase, there is 
a unique, singlet ground state with an energy gap to the first exited state. For = 6, 
this is also what we have here, and hence one could expect that may be similar to 
ip^ in this parameter regime. In fact, the overlap between the two states is seen to be 
very close to unity in an interval within the Haldane phase. The maximum is 0.9990 
and occurs for = —0.3213. A nonzero overlap is also seen in the dimer and trimer 
phases, where the condition ^ . Sj = is fulfilled. In the ferromagnetic phase, on the 
other hand, the overlap vanishes because the ground states of (I155p all have total spin 
A^. The results presented below for the Renyi entropy and the two-point spin correlator 
suggest that the spin 1 model at level = 2 is gapless for A^ — )■ oo, and hence we do 
not expect the two models to be close to each other in that limit. 



N 




(155) 



1=1 
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4.3. Entanglement entropy and two-point spin correlation function for j = 1 

It is interesting to ask, whether there is a simple connection between the properties of 
the conformal fields used to construct the wave function and the properties of the wave 
function. This turns out to be the case for the uniform HS model, where it is known 
that the low energy states are described by the SU{2)i WZW model in the large N limit 
[3] . To investigate whether a similar result could be true for the spin 1 case, we compute 
below the Renyi entropy and the two-point correlator of the spin 1 chain numerically 
and compare the results to those expected for the SU{2)2 WZW model. 
We start with the Renyi entropy 

^-ln[Tr(pi)], (156) 

where pi is the reduced density operator of the first L spins obtained by tracing out 
the state of the last N — L spins. Rewriting this equation, we can express the entropy 
in terms of an expectation value 

XI XI lV'«i,...,«^(2:i,---,^7v)lV4v..,^'^(^b--->2;7v)|^ (157) 



o{2) 



X 



X \'^s^,...,sAzU---.Zn)\' 



-2 



which we evaluate numerically by use of the Metropolis Monte Carlo algorithm [6] . The 
Renyi entropy of a conformal field theory with central charge c is [291 EOj [31] 

5f ) = ^ ln[sin(7rL/A^)] + constant, (158) 

and c = 3/2 for the SU{2)2 WZW model (recall (II])). In figure El we have plotted Sf^ 
as a function of ln[sin(7rL/iV)]/4, and indeed the points fall almost on a straight line. 
Considering c as a fitting parameter, we find that c = 1.395 is the best choice to describe 
the data, but c = 3/2 also gives a fit that approximately follow the data points. 

Let us next consider the two-point spin correlation function. For the spin 1 model at 
level = 2, we do not obtain a closed set of equations for the two-point spin correlators 
from (13^ and fl6Up . so we use Monte Carlo simulations instead. Results for the uniform 
model and = 200, 500 and 1200 are shown in figure [6] Affleck found that the spin 
correlation function for an infinite chain of spin 1 fields in the WZW model at level k = 2 
is proportional to (— l)"n~^/^ for large n [32], possibly with a multiplicative logarithmic 
correction of the form (A + i?ln(n))^/^, where A and B are constants [2Q\- Referring 
back to f ll33p and also to the scaling behaviour of the correlator in the dimer model 
discussed in section 13.3.41 it is very natural to expect that the most significant term of 
the correlation function for the finite spin 1 chain is proportional to (— l)"[sin('7rn/A^)]^ 
in the asymptotic region n ^ 1, where the exponent b may or may not be equal to —3/4. 
We therefore plot ln[(— l)"(tQt")] versus ln[sin(7rn/A^)] in figure El This transformation 
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Figure 5. Renyi entropy of the state (|148l) for zj = exTp{2TTij / N) and N = 200 
obtained from Monte Carlo simulations of exp(— 5')^ ). The data shown is the average 
of 40 trajectories with different initial conditions, and the variation of the mean of 
these trajectories is used to estimate the error bars. The solid line is a fit of the form 
S^j^^ = cln(sin(7ri/iV))/4 -|- a, where a and c are fitting parameters, and the dotted 
line is a fit of the same form but with c fixed at 3/2. 
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Figure 6. The functional dependence of the two-point spin correlation function on 
separation of the spins for Zj = exp(27rij/A^) and N — 200, 500 and 1200. The 
points are obtained from Monte Carlo simulations, the solid lines are fits of the form 
ln((— = a + feln[sin(7m/A^)], where a and b are fitting parameters, and the 
dotted lines are fits of the same form but with b fixed at —3/4. As in the previous 
figure, the error bars are estimated from the variation of the mean of trajectories with 
different initial conditions. 
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indeed gives data lying almost on straight lines despite the fact that we are considering 
chains of finite length. Using a fit to extract the slope of the lines, we find that the 
exponent is around —0.69. A similar exponent was found for the state considered in 
[33] . As in [33], we note that the deviation of b from —0.75 may be due to the presence 
of a logarithmic correction since a logarithmic correction of the type given above is 
equivalent to a power law to first order in the exponent: ra^/^ (1 + In(n)e)^/^ for e 
small. Altogether, it is hence possible that the properties of the wave function are those 
of the SU {2)2 WZW model, but it may also be that there are some deviations. 

5. Conclusion 

In conclusion, we have used a specific set of null vectors to derive quantum spin 
Hamiltonians, whose ground states are the chiral correlators of the SU {2)k WZW model, 
and we have shown that the same null vectors lead to a set of algebraic, linear equations 
relating the spin correlators within each model. At level /c = 1, these equations are 
sufficient to compute both two-point and higher-point spin correlators simply by solving 
a set of linear equations, which can even be done analytically in some cases. The models 
involve Hamiltonians with long-range interactions, and explicit analytical expressions 
for the corresponding ground state wave functions can be found in many cases. The 
examples presented in the paper include both new spin models and new results for 
already known models. 

It is interesting to note that our construction in some sense resembles that of 
the AKLT [21 |37] bond in matrix product states (MPS) and frustration free parent 
Hamiltonians. As in MPS, the components of the wave functions are constructed 
from a product of N operators (see ( l34l) ). where the jth operator depends only on 
the state of the jth spin, but the operators are chiral conformal fields instead of finite 
dimensional matrices. Similarly, the Hamiltonian is a sum of Hamiltonians (see (|53|) ) 
that annihilate the ground state individually, but the Hamiltonians involve long-range 
interactions rather than only local interactions. These differences, in fact, circumvent 
one of the limitations of MPS, namely how correlation functions may decay algebraically, 
and thus one can also use the states to describe critical behaviour. 

The uniform Haldane-Shastry model is closely related to the Calogero-Sutherland 
(CS) model, whose Hamiltonian contains a kinetic term plus an inverse square coupling 
between particles moving on a circle [M] (see [10] for a review). In turn, particular 
eigenfunctions of the CS model are in direct correspondence with the correlation 
functions in boundary conformal field theory. In fact, the CS Hamiltonian can be 
obtained from the null vector decoupling equations at level 2 of boundary fields and a 
bulk primary field at the centre of the disk [35l [36] • These results suggest a possible 
relation of our work to those of [351 ES] , although the null vectors used there and those 
used in the present work are different. 

The present paper raises the question, whether the properties of the wave functions 
are directly related to the properties of the fields used to construct the wave functions. 
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The properties of the uniform model with k = 1 and j = 1/2 [3] and the numerical 
results presented above for the two-point spin correlator and the Renyi entropy for the 
uniform model with k = 2 and j = 1 suggest that the properties of the wave function 
may be similar to, but not necessarily identical to, the properties of the fields. We are 
currently investigating this question in more detail. 

There are several ways in which the explicit examples given in the paper can be 
generalized. The construction in section |2] already includes the possibility to consider 
models at higher levels, models in which the fields used to construct the wave functions 
do not all have the same spin index, and models in which the individual terms H^^^ j are 
given different weights in the final Hamiltonian. Furthermore, we have only considered 
one particular set of null vectors, but other Hamiltonians can be derived by starting 
from other null vectors. The present work shows that translationally invariant models 
of one- dimensional spin chains with long-range interactions can be achieved by choosing 
the 2;j-coordinates to be uniformly distributed on the unit circle. More generally, the Zj- 
coordinates are numbers in the complex plane, and various interaction patterns can be 
achieved by adjusting these coordinates. This observation facilitates the construction of 
higher dimensional spin models. In particular, it is natural to describe a two-dimensional 
spin lattice in the above framework by interpreting, e.g., the real and imaginary parts 
of Zj as the coordinates of the physical position of the jth spin in the lattice. This 
construction will be presented elsewhere [3H]- Finally, the idea of finding Hamiltonians 
from null vectors can also be applied to any Lie group and to WZW models based on 
supergroups. Altogether, the present work thus opens up possibilities to find a huge 
class of spin models. We conclude that null vectors are a powerful tool to gain insight 
into the behaviour of many-body systems, and we believe that more results will be 
obtained along these lines. 

Note added in proof. A spin 1 model with the same type of interactions as the model discussed in 
section 4 has been derived for the uniform case in j39j , and an application of conformal field theory to 
propose wave functions for the fractional quantum Hall effect is considered in |40) . 
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Appendix A. Alternative derivation of the generalized HS model 

In [6], it was found by use of the Knizhnik-Zamolodchikov equation that the state fl67|) 
is an eigenstate of the generalized HS Hamiltonian fl87|) . We provide the details of this 
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derivation in the following. We start from the KZ equation 

|^^=E^^' ^ = l,---,^, (A.1) 

for the chiral correlator ( l62l) with components ( |67l) . Using the conformal transformation 
z = e'^, one can map the complex plane into the cylinder, where the chiral correlator 
becomes 

1pcyl{Wl, Wn) = (^^~^ ■ ■ ■ {^^^ ^plane(^l, ■ ■ ■ , Zn) (A. 2) 

and /i = 1/4 is the conformal weight of the spin 1/2 primary field. Expressing ipcyi in 
terms of the Zi variables and using ( 1671) for ipp\ane{zi, . . . , zn) = ip{zi, . . . , z^), one has 

1pcyl{Zl, . . . , Zn) = Y\_ ^plane(2;i, ■ ■ ■ , Zn) = Y\^^^ ~ ZjY'"'^'^ Y\_ ^i^' (^--3) 

i i<j i 

The KZ equation for this wave function is 
d 

where Wij = {zi + Zj) / {zi usually. 

Using the explicit expression for the wave function in ( lA.SI) and the relations 

J2 = = -1, = 0' (A-^) 

one finds 

= U E ^/^^^ + i = \ E - ^1^) ....^c.i.(A.6) 

Simplifying the above expression one gets an abelian version of the KZ equation 



Taking another derivative, one finds 

(d \^ 1 . z-z- 1 . 

^ig^, j ^cyi = -2^2^ 16 5Z ^^J'^^fc ("^-^^ 

where we have used = 1 and 

Zi^Wij = j^. (A.9) 

dz, 4 

Summing over i in the last expression, one gets 



'\Y^ + 1^ 5Z 5Z 5Z ^^J^^^ ^^'^^ ^=yi- 
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The triple sum can be written as 

so that 

^^cyi = - 2 1] ^ ^i^j ^=yi + 16 + 16 5Z ^^i^^fc ^^'^^ ^^yi- ^^-^2) 

Exchanging i ■v^ k, this becomes 

^^cyi = "2 5Z ^ ^^^^^ + 16 5Z ^'^^ ^ 16 ^ ^ ^^'^^^ 
Using (EZD, ([79]) and s^Sj = -A^ in f OTTSjl . we find 

^ = E ^ + 1^ E + 1^ E ^^^^-[^ - 2 + + -.-(Q - c,)] (A.14) 



1 / ZiZj 1 2 \ 1 2 N — 2 ^r~^ 

= - 2 E (^1^ - i-.- J + E + ^ E 
+^ E -'^^^ - " ^ E -'^-'^ + ^^1^ E ^^^j- + ^ E 
+1^ E - 9) = + E + 1^ E - c,)- 

From (lA.7p . we have 

1 V d 1 V 

- ^ ^cyl = Ye E ■5i'SjWij(Ci - Cj)lpcyU (A.15) 

i * i^j 

SO finally 

- ^ E ^^^^1") ^^yi = (-1^ + 1^ E ^j) ^cyi. (A.16) 

The operator on the left hand side of f lA.16P can also be computed from the KZ 
equation (1A.4I) . Taking the derivative of (1A.4I) and following the same steps as before, 
one finds 

9 AV^eyi = -6 E ^t^tl^^y^ + E wlm)'^^y^ + E E «^H«;.,t^t?tti,Vcyi. (A.17) 



The following identities are needed 

3 1 
16 ~ 2 



mf = T^-ottt';^ 1^3. (A.18) 



= {t-f; + ea,ctlt\t% t^j^k, (A. 19) 

tXj = T^T" - ^N, (A.20) 
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which finally lead to 

2 ZiZi „ „ 1 9 A^(A^ - 2) 1 , , „ „ 



3 ^ 



where we have used T"'T"'i/jcy\ = because ipcyi is a singlet. It also follows from flA.4p 
that 



1 x d 1 V 



i 

and hence 

^-^E^^-^l;)^' 



(A.21) 



cyl 



IE 



Ar(A^-2) 
48 



V'cyl- 



Equating flA.16P and flA.22p . we arrive at 



E 



-E' 



Ar(Ar + 1) 

16 ^ 



16 



cyl- 



(A.23) 



Since ipcyi only differs from ip by an Si-independent factor, we can replace V'cyi by ip in 
f lA.23p . and it follows that is an eigenstate of ( 187|) as claimed. In [6], it was found 
numerically for small systems that ip is, in fact, the ground state of (187|) . With the 
results of the present paper, however, it is clear from the construction that ip is the 
ground state. We note also that the above derivation is not immediately generalizable 
to k> I. 



Appendix A.l. The Inozemtsev invariant 

Let us finally use the above approach to derive the generalization (19T]) to the nonuniform 
case of the Inozemtsev invariant. We start from the second derivative 

(d \^ 1 . z-z- 1 . 

^i-Q^] ^cyi = "2 E ^ ^'^^ ^'^^+16 5^ 5Z ^^J'^^fc ^cyl(A.24) 

and take another derivative 



^cyl V'cyl \ dZi) 2^, dZi\ zfj 

E E 1? ^^^^i^fc + 3^ E E E (^^i^ifc) s.sfc 

+ ^ E E E E ^ij^ikWil SiSjSkSi. 

The terms involving 4 different spin variables vanish by virtue of the equation 

WijWikWii + (permutations of k,l) = 0, i ^ j ^ k ^ I. (A. 26) 
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Then, using the following identities 

d 1 

Zk-^{wkiWkj) = -{wki - Wkiwlj + Wkj - Wkjwli), (A.28) 

we find after several manipulations 

3 V 

AsV^cyl = 5Z ^ij^i^ji^cyh Ci,2 = (A.29) 

Hence, 

(^3 + ^E^^-2^^|-jv^cyl = 0. (A.30) 

The left hand side can be computed using the KZ equation. After a long computation 
one gets H^ipcyi = and hence H^ip = 0, where is given by fPTl) . 

Appendix B. Hamiltonian for k = 2 and j = 1/2 

The Hamiltonian at level k = 2, which has the chiral correlator (0i/2,si(-2i) • • • 4>i/2,SMi^N)) 
of spin 1/2 fields as its ground state, follows from f l32|) . fj45l) . fl50l) and fl53l) . Assuming 
l^il = 1 Vz, we find 

+ 2 E ^lMi2Vi2 + l E Wii.wu.Wii.Wii.t'l^tltlt]^ 

2 ^ 

+ 5 E Wii,Wii2?i;ii3(4wii, - Wii^ - Wii.^t1t1t\t\ 

and 



1 \ 2 r 

+ 2^^1*3^1*3 l^^l ^^2 + 5 E E ^i^in^i^U2W^«3^ii4 



n7^*2^«37^M j(7^n,j2,j3,M) 



j.a j.a J.6 J.6 
''n''i2''i3'^«4- 



The spectrum of (\B.2\i for the uniform case and A^ = 6 is given in table IB1[ Since 
two 01/2 fields fuse to a 0o and a 0i field, and since the fields should altogether fuse 
to 00 to give a nonzero correlator, the ground state is 2^/^-1 times degenerate. Apart 
from this degeneracy, all degeneracies between spin multiplets can be explained from 
the invariance of the chain under translation by one site and under a fiip of the chain. 
This suggests that there is no Yangian symmetry at level k = 2. Another difference 
compared to the uniform HS model is that not all energies are integers. 
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Table Bl. Spectrum of the spin 1/2 Hamiltonian (|B.2p for zj = exp(27rij/A^) and 
A'^ = 6 sites. E is the energy, S is the spin, p„ is the momentum, and ^ = 25 + 1. 



E 


S 


3p„ 

TT 


# 


table continued 






table continued 






0.0000 





-1 


1 


20.0000 


1 


-1 


3 


107.2000 


3 





7 










1 




1 


1 


3 


111.2000 


2 - 


1 


5 







1 


1 


34.4000 


2 


-2 


5 




2 


1 


5 







3 


1 




2 


2 


5 


117.1314 


1 - 


2 


3 


0.2715 


1 





3 


47.2000 


2 


3 


5 




1 


2 


3 


2.8686 


1 


-2 


3 


52.0000 


1 


3 


3 


160.0000 





3 


1 




1 


2 


3 


58.9285 


1 





3 
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